Quadratic interaction effect on the dark energy density in the universe
Deveci Derya G1, 2, Aydiner Ekrem1, †
Department of Physics, Faculty of Science, İstanbul University, Fatih Tr-34134, İstanbul, Turkey
Department of Opticianry, Altlnbaş University, Tr-34144 İstanbul, Turkey

 

† Corresponding author. E-mail: ekrem.aydiner@istanbul.edu.tr

Abstract

In this study, we deal with the holographic model of interacting dark components of dark energy and dark matter quadratic case of the equation of state parameter (EoS). The effective equations of states for the interacting holographic energy density are derived and the results are analyzed and compared with the solution of the linear form in the literature. The result of our work shows that the value of interaction term between dark components affects the fixed points at far future in the DE-dominated universe in the case of quadratic EoS parameter; it is a different result from the linear case in the theoretical results in the literature, and as the Quintom scenario the equations of state had coincidence at the cosmological constant boundary of –1 from above to below.

1. Introduction

Recent Supernova Type Ia observations (SNIa)[1,2] have showed that our Universe appears to be accelerating and the universe was expanding more slowly than it is today. They also suggest that there is a form of energy which is unknown and is shelving to accelerated expansion of the universe by energy contribution ΩDE ≃ 0.75. This draws attention in the astronomy and particle physics research. Besides, cosmic microwave background (CMB) observations[3] provide evidence of the Universe to be expanding at an increasing rate. Cosmological observations suggest that a large part of the Universe arises from dark components which are named dark energy (DE) and dark matter (DM).[46] There are a number of dark energy candidates such as the cosmological constant model and scalar fields model including quintessence scenario, k-essence model, Chaplygin gas, holographic dark energy (HDE) models, and phantom dark energy model.[7] Also two famous cosmological constant problems survive that researchers try to solve: the fine tuning and coincidence problems. In the observational cosmology, the interaction between dark components giving EoS parameter of DE ωDE = pDE/ρDE plays an important role in solving the mentioned problems where pDE is the pressure and ρDE is the energy density. According to literature, in order to derive accelerated expansion of the Universe the EoS parameter of DE must satisfy ωDE < −0.33 condition.[8] The EoS of the positive cosmological constant represented by Λ is the simplest candidate of the DE which has value −1, the value for quintessence model is −1 < ωDE < −0.33, for k-essence model is ωDE < −1 or ωDE > −1. In the Chaplygin gas and HDE model ωDE > −1, and not crossing −1. In the phantom case ωDE is smaller than −1. Most cosmological models assume that dark components only interact gravitationally.[9] In the recent works the parameter of EoS evolves across −1, the cosmological constant boundary of this phenomenon is called the quintom scenario, that was first pointed out in a pioneer work[10] and was later comprehensively reviewed in Ref. [11]. In cosmological models the idea that DM and DE are not evolving separately but interact with each other was introduced first to justify the currently small value of the cosmological constant but they were found to be very useful to alleviate the problem of coincidence namely: Why are DE and DM densities of precisely the same order today? discussed in recent works.[1214] Since DE and DM dominate the energy content of the universe today, it is equally reasonable to assume that these dark components could interact among themselves and with other components. The mutual interactions between DE and DM can be suggested as a possible scenario of the Universe evolution. On the other hand, the interacting dark components are provided a new direction to understand the nature of DE. The particle physics can be explained by an effective field theory with a (short distance) UV cut-off that is less than the Planck mass given by Mp. For an effective quantum field theory in a box of size L with UV cut-off Λ the entropy S scales extensively, SL3Λ3. According to Bekenstein[15] the peculiar thermodynamics of black holes, which have the maximum entropy in a box of volume is L3, and they behave non-extensively, therefore Bekenstein entropy bound is called Bekenstein–Hawking bound. The recent works[16] show that where SBH caLled by the entropy of a black hole with radius L acting as an IR cutoff, however it cannot be chosen independently of the UV cut-off, and scales as Λ3. It means that UV cut-off is connected to the IR cut-off (long distance) and the density of zero point energy is proportional to UV cut-off Λ. The whole system total energy having size L may not exceed the mass of black hole with the same size. It means that the black hole entropy is directly related to the area of the event horizon. The holographic principle[17] was originated by black hole thermodynamics. It states how quantum theory and gravity interact to construct the reality that we are in and it suggests that the Universe can be seen as two-dimensional (2D) information on the event horizon and can be thought of as encoded on a lower-dimensional boundary to the region, therefore IR cut-off related to the DE is the size of the event horizon. HDE[18] arises from a theoretical attempt of applying the holographic principle to the DE problem; it has been introduced recently.[19] In this model the future event horizon is chosen as the characteristic length scale L and the DE density given by , where c is a free dimensionless parameter. Many research studies try to solve the mentioned cosmological problems with working different interacting models. Recent observations of clusters of galaxies indicated that the density of matter is less than critical density. Observations, especially CMB signs, specify that the Universe has flat form and the Universe is dominated by DE with negative pressure driving to accelerate expansion of the Universe. Research where the interacting DE models are derived with linear form of EoS parameter, concluded that the mutual interactions of DE and DM may play a crucial role in the dynamics of the Universe by alleviating the mentioned cosmological problems.[2024]

2. Interacting model

Generally it is assumed that interaction between DM and DE can be defined as

with the total energy density ρtot = ρDE + ρDM where ρDM and ρDE are energy densities of DM and DE respectively. pDM and pDM are also the pressure of DM and DE, respectively. Q represents the overall conservation of the energy–momentum tensor indicating a transfer from DE component to DM component and vice versa. In the literature there are non-zero matter components of ρDM[24] but in this study pressure is taken as zero for the matter component (pDM = 0) because we know that the normal dust matter is pressureless. In this study pDM is considered as a specific case of a general quadratic EoS of the form[25]
We investigate a quadratic EoS parameter focusing on the quadratic case and ignore higher order terms. In high-energy regime p0 can be chosen as zero. The parameters A1 and A2 set the characteristic energy scale. They can be chosen to be the same. For the standard linear EoS parameter for a cosmological fluid, p0 and A2 set to be zero, but now we are interested in the A1 = A2 condition with ωDE value. Rewrite Eq. (2) as
This is decaying of the HDE into cold dark matter (CDM). We assume a form of interaction Q = ΓρDE with the decay rate Γ. Equation (1) can be reorganized by using Eq. (3)
where the effective equations of states are defined as follows:
in this equation the ratio of the two energy densities is given by u = ρDM/ρDE. The homogenous and isotropic universe can be easily described by Friedmann–Robertson–Walker (FRW) metric[26]
where k refers to the curvature of the space with k = −1, 0,1, the curvature of space depending on whether the shape of the Universe is an open 3-hyperboloid, flat (i.e., Euclidean space) or a closed 3-sphere, respectively, and a is the cosmic scale factor depending on time and is a measure of the size of the Universe representing relative expansion of the Universe. x = 0 when a = 1 since x = ln(a). In order to combine the energy density with the curvature of the Universe, we employ the first Friedman equation given by
where ρtot = ρDE + ρDM, c is velocity of the light, c = 1 is known for the flat universe. is the reduced Planck mass and H = ȧ/a is the Hubble parameter. The geometry of the Universe is defined by the density parameter that is expressed in terms of the ratio of the Universe density to the critical density. This critical density is also called closure density which is a remarkably small number for our Universe and it is described by ρcr. Ω0 is defined as the average matter Universe density divided by a critical value of that density (Ω0 = ρ0/ρcr). This selects also one of the three possible geometries depending on whether Ω0 is equal to, less than, or greater than 1 and it governs whether the curvature is negative (Ω0 < 1), positive (Ω0 > 1), flat (Ω0 = 1). Today, the research shows that the Universe is nearly flat, therefore the average density is nearly close to the critical density (ρcr = 10−26 kg/m3). Define the dimensionless density parameters for DE, DM, and the curvature relative to critical density parameter;
where is the critical energy density and the First–Friedman equation can be rewritten by using Eqs. (7) and (8)
which is for the holographic condition involving the future event horizon and then expresses the density ratio in terms of Eq. (8) by using Eq. (9)
where u = ρDM/ρDE which should be a solution for the coincidence problem. Its derivative can be rewritten by using Eq. (10)
which means that the ratio of the densities depends on the defined form of the interaction. In the evolution when is equal to 0, which means that the ratio of the energy densities provides constant value for the DE and DM. means that the dark energy dominates the system and it is not vanishing. If the vanishing dark energy density is being apparent. The IR cut-off of the Universe in the case of the holographic model with the HDE density is defined in terms of L
where L is IR cut-off which can be identified by the radius of the future event horizon r defined as follows:
while referring holographic model, when DE-dominated universe model is assumed, the future event horizon has a tendency to a constant value with the order H−1, therefore equation (13) can be rewritten by L = 1/H. In this study we have chosen the radius of r(t) as follows:
where the event horizon is given by Rh and the radial coordinate is given by
using these definitions ΩDE and ρcr, we can easily obtain
using Eqs. (14)–(17) we gain
and then using the definition Ωk and Eq. (18) we obtain
where cos y = cos y for k = 1, cos y = y for k = 0, and cos y = cosh y for k = −1 with . If equation, (13) is used cos y can be rewritten as . Taking again the holographic energy density definition ρDE by using Eq. (4) one finds the holographic energy equation of the state by using Eq. (13) takes the following form
where . In this work, the decay rate can be chosen by Γ = 3b2(1+u)nH (n is a constant) where b is the coupling constant and 1 ≥ n. Substituting this relation into Eq. (5) with the chosen decay rate we find the effective state equation
Now we derive here two dimensionless coupled equations with solutions determining the effective equations of state. ΩDE is coming from the derivative of u given by ΩDM/ΩDE and Ωk is coming from the derivative of the rk given by Ωk/ΩDE. The corresponding coupled equations are derived by using Eqs. (11), (12), and (21)
where x = ln (a/a0) and . When we compare our results with literature it can be seen that our model is in accordance with the literature.

3. Numerical results

In this study we have obtained the state equation for the interacting HDE density by using quadratic form of the EoS parameter. The numerical results are analyzed for the cases of k = 0, 1,−1 the curvature of space flat, closed, and open forms of the universe, respectively. In this study we have used the effective equations of state (, ). To analyze the results of our numerical solutions we represent the behavior of the energy densities and effective equations by plotting versus x values in the case of quadratic form of the EoS parameter.

In the literature the effective equation of the state for n > 1 diverges for small ΩDE for choices of the holographic condition that could be consistent with observations, therefore in literature the parameters were assigned as c = 1, n = 1, and b2 = 0.1 and 0.2 with the initial conditions at x = 0, ωDE(0) = 0.72, ωk(0) = 0 for k = 0, ωk(0) = 0.08 for k = 1, ωk(0) = −0.01 for k = −1.

First of all the evolution behavior of the effective equation of state (blue) and (red) in Eqs. (21) and the evolution of ΩDE and Ωk with respect to x is qualitatively determined in Eqs. (22) are plotted in Figs. 1, 2, and 3 to analyze the numerical results and compare them with previous literature.

Fig. 1. (color online) (a) The effective equations of state (blue) and (red) versus x where x = ln(a/a0) and (b) the density parameters of ΩDE (blue) and Ωk (red) versus x are represented for the non-flat universe (closed case) b2 = 0.2, for an interaction of n = 1 with the initial conditions ΩDE(0) = 0.72, Ωk(0) = 0.08 where k = 1.
Fig. 2. (color online) (a) The effective equations of state (blue) and (red) versus x where x = ln(a/a0) and (b) the density parameters of ΩDE (blue) and Ωk (red) versus x are represented for the flat universe b2 = 0.2, for an interaction of n = 1 with the initial conditions ΩDE(0) = 0.72, Ωk(0) = 0.0 where k=0.
Fig. 3. (color online) (a) The effective equations of state (blue) and (red) versus x where x = ln(a/a0) and (b) the density parameters of ΩDE (blue) and Ωk (red) versus x are represented for the non-flat universe (open case) b2 = 0.2, for an interaction of n = 1 with the initial conditions ΩDE(0) = 0.72, Ωsk(0) = 0.01 where k = −1.

In Fig. 1 the density parameter and effective equations are plotted with the initial parameters b2 = 0.2, for an interaction of n = 1 with the initial conditions ΩDE(0) = 0.72, Ωk(0) = 0.08, where the positive curvature constant k = 1. In this case the density is so high that the gravitational attraction can stop the expansion and collapse backward named “big crunch”. This kind of universe is described as being a closed universe, or a gravitationally bound universe with a small positive curvature (Ωk = 0.02). As shown in Fig. 1(a) the two effective equations of state start different points. In far future two effective equations take the same value () between 0 and −1, which means that there is a stable point in the far future. It is clear in Eq. (12) when r is equal to zero ( meaning that the Universe may be re-collapsed in the future time. From Eq. (9) the density of matter exceeds the critical value (0.28) when the curvature term goes to zero at the far future it is in accordance with Ref. [20]. In Fig. 1(b) it is seen that Ωk has maximum positive curvature with value 0.04 and Ωk goes to zero for the past and far future and the density parameter of ΩDE is a monotonically increasing function of x which goes to 0.77 at far future, therefore ΩDM = 0,27. The research shows that the observed best fit of the closed model of the Universe has ΩDM = 0.415, ΩDE = 0.630.[30]

In Fig. 2 the density parameter and effective equations are plotted with parameters b2 = 0.2, for an interaction of n = 1 with the initial conditions ΩDE(0) = 0.72, Ωk(0) = 0.0 where the curvature constant k = 0. In theory a space with zero curvature is called a flat space which is non-compact and extends infinitely far in any direction and in this case the density value is equal nearly to the critical value, so the Universe may expand forever. As shown in Fig. 2(a) the two effective equations of state start also from different points. In the far future the two effective equations approach their limiting (equal) value () which means that there is a stable point in the far future. The two effective equations coincide at −0.98 fixed point. In Fig. 2(b) Ωk has zero value for the past and far future. The density parameter ΩDE has monotonically increasing behavior like non-flat universe 0.78 at the far future; it is compatible with observations.[27]

In Fig. 3 the density parameter and effective equations are also plotted with the same parameters except the initial condition Ωk(0) = −0.01 where the Universe has open form with negative curvature, space has infinite volume. This space–time also has space expanding forever in time. In Fig. 3(a) in the far future the two effective equations are also coincident at a fixed point (), and take the value between 0 and −1. In Fig. 3(b) Ωk has negative value with tending zero value for the past and far future. The recent research shows that, in the case a non-flat universe, the models with small negative value of Ωk are a better fit than ΛCDM model.[30] The density parameter ΩDE is also monotonically increasing.

In Fig. 4 the equation of states and density parameters are plotted with different interaction parameters to see the effect on the nature of the holographic interaction. In the case of b2 = 0 non-interaction we see ΩDE(x) → 1 as x value increases. When we have chosen the interaction parameter value to be 0.1 as shown in Fig. 4(a) in the far future two effective equations are also coincident at a fixed point, and take the value between 0 and −1, however it is a bigger value than the b2 = 0.2 case. In Fig. 4(b) ΩDE(x) is also a monotonically increasing function of x and reaches 0.90 value at the far future.

Fig. 4. (color online) (a) The evolution behavior of the effective equations of state (blue) and (red) versus x where x = ln(a/a0) and (b) the evolution behaviors of the density parameters of ΩDE (blue) and Ωk (red) versus x are represented for the non-flat universe (open case) in the case of quadratic EoS with the parameters b2 = 0.1, for an interaction of n = 1 with the initial conditions ΩDE(0) = 0.72, Ωk(0) = 0.08 for k = 1.
4. Conclusions

In cosmology, driving to an accelerated universe, a quadratic EoS parameter has been shown to be concerned with the expounding of DE and DM,[28,29] it is also appearing in the energy density evolution equations and may be liable for the annihilation of anisotropy and inhomogeneity. In contrast, for the standard general relativistic cosmology case, a linear EoS parameter represents the high isotropy of the Universe. In this study, we have chosen the DE-dominated non-flat universe and three parameters b2 (interaction strength), c, n are taken in the observational ranges ΩDE, ΩDM, Ωk. First of all we have extended the interacting HDE model to the non-flat model of the Universe with quadratic EoS parameter which represents the state of the composition of two comparable components that arise from decaying of HDE into CDM. Using HDE, we have obtained the equation of state for interacting HDE density enveloped by L universe. The density parameters of ΩDE and Ωk have calculated numerically for the quadratic EoS parameter, and then the effective equations and are calculated. The behavior of the density parameters and effective equations are analyzed. As can be seen from the numerical results mentioned above the effect of the non-flat universe is trivial, because from Figs.1, 2, and 3(a) it can be seen easily that Ωk → 0 for the past and future, and this means that for the non-flat universe with value k ≠ 0, the phantom regime could not be caused for which , even if one comprises an interaction between HDE and DM where . For the non-interacting case of b2 = 0 we see that ΩDE(x) → 1 as x increases. As shown in Figs.1, 2, 3, and 4(b), ΩDE(x) is a droningly increasing function of x with the interacting cases of b2 = 0.1 and b2 = 0.2. The DM density is affected by small changes of the interaction term, but the DE density is not. This means that the interaction term affects strongly the effect on the nature of the holographic interaction in the DE-dominated universe in the case of the quadratic EoS parameter, which is different from the linear case. These figures show also that the quadratic effective EoS of the two equation states of and will take negative values which are greater than −1 during the whole evolution of the Universe. This means that the holographic energy density preserves its nature for the future even if the interaction term is changed and as the Quintom scenario the equations of state had coincidence at the cosmological constant boundary of −1 from above to below. As can be seen in Fig. 4 a small change is important when it comes to quadratic additions to EoS parameter. This refers that the interaction term affects the fixed points at the far future in the DE-dominated non-flat universe. The effective equation of matter component is also affected by small changes in the interaction parameter; however the effective equation of DE is not affected strongly. When it is compared with the recent research, it is concluded that use of quadratic EoS parameter leads to the overriding of the term of the interactions between DM and DE, therefore the effective equation of state parameter fails to evolve from the regime greater than −1 to the regime less than −1 or the opposite way. In this regard, this can lead to the conflict between the theoretical and observational expectations.

Reference
[1] Riess A G Filippenko A V Challis P et al. 1998 Astron. J. 116 1009
[2] Perlmutter S Aldering G Goldhaber G et al. 1999 Astron. J. 517 565
[3] Fabris J C Goncalves S V B De Souza P E 2002 arXiv: preprint astro-ph/0207430
[4] Setare M R 2007 J. Cosm. Astron. Phys. 2007 023
[5] Jussi V et al. 2015 J. Cosm. Astron. Phys. 2 015
[6] Yang W Li et al. 2016 J. Cosm. Astron. Phys. 2016 007
[7] Copeland E J Mohammad S Tsujikawa S 2006 Int. J. Mod. Phys. D 15 1753
[8] Sadeghi J Setare M R Banijamali A Milani F 2008 Phys. Rev. B 66 92
[9] Caldera-Cabral G Maartens R Urena-Lopez L A 2009 Phys. Rev. D 79 063518
[10] Feng B Xiulian W Xinmin Z 2005 Phys. Rev. B 79 063518
[11] Cai Y F Saridakis E N Setare M R Xia J Q 2010 Phys. Rep. 493 1
[12] Hu B Ling Y 2006 Phys. Rev. D 73 123510
[13] Zimdahl W Y 2005 Int. J. Mod. Phys. D 14 2319
[14] Dalal N Abazajian K Jenkins E Manohar A V 2001 Phys. Rev. Lett. 87 141302
[15] Bekenstein J D 1994 Phys. Rev. D 49 1912
[16] Cohen A G Kaplan D B Nelson A E 1999 Phys. Rev. Lett. 82 4971
[17] Bousso R 2002 Rev. Mod. Phys. 74 825
[18] Wang S Wang Y Li M 2016 arXiv: 1612.00345
[19] Wang B Gong Y Abdalla E 2005 Phys. Rev. B 624 141
[20] Kim K H Lee H W Myung Y S 2007 Phys. Rev. B 648 107
[21] Kim H Lee H W Myung Y S 2006 Phys. Rev. B 632 605
[22] Wang B Lin C Y Abdalla E 2006 Phys. Rev. B 637 357
[23] Wang B Abdalla E Atrio-Barandela F Pavon D 2016 arXiv: 1603.08299
[24] Berger M S Shojaei H 2006 Phys. Rev. D 73 083528
[25] Aydiner E 2016 arXiv: 1610.07338
[26] Setare M R 2006 Phys. Rev. B 642 1
[27] Sheykhi A 2009 Phys. Rev. B 680 113
[28] Kishore N A Marco B 2006 Phys. Rev. D 74 023524
[29] Sharma R Ratanpal B S 2013 Int. J. Mod. Phys. D 22 1350074
[30] Spergel D N Bean R Doré O et al. 2007 Astron. J. Suppl. Ser. 170 377